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In our recent paper entitled "Pairing mechanism of high-temperature superconductivity: Experi- 
mental constraints" (to be published in Physica Scripta), we review some crucial experiments that 
place strong constraints on the microscopic pairing mechanism of high-temperature superconductiv- 
ity in cuprates. In particular, we show that phonons rather than spin-fluctuation play a predominant 
role in the microscopic pairing mechanism. We further show that the intrinsic pairing symmetry 
in the bulk is not d-wave, but extended s-wave (having eight line nodes) in hole-doped cuprates 
and nodeless s-wave in electron-doped cuprates. In contrast, the author of the Comment (to be 
published in Physica Scripta) argues that our conclusions are unconvincing and even misleading. In 
response to the criticisms in the Comment, we further show that our conclusions are well supported 
by experiments and his criticisms are lack of scientific ground. 
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In our recent paper [![, we review some crucial ex- 
periments that place strong constraints on the micro- 
scopic pairing mechanism of high-temperature super- 
conductivity in cuprates. In particular, we show that 
phonons rather than spin-fluctuation play a predominant 
role in the microscopic pairing mechanism. We further 
show that the intrinsic pairing symmetry in the bulk of 
cuprates is not d-wave, but extended s-wave (having eight 
line nodes) in hole-doped cuprates and nodeless s-wave in 
electron-doped cuprates. However, Plakida 2] has raised 
strong criticisms on these conclusions based on an over- 
simplified polaronic model and some experimental results 
that have been misinterpreted. Below we will show that 
our conclusions are well supported by experiments and 
the criticisms raised in the Comment are lack of scien- 
tific ground. 

In the Comment 0], the author first considers the 
oxygen-isotope effect (OIE) on T c by taking into ac- 
count the observed oxygen-isotope effect on the in-plane 
effective supercarrier mass. He has used the weak- 
coupling BCS-like formula [Eq. (2) in the Comment] to 
calculate the doping dependence of the OIE on T c in 
La2- 2: Sr :r Cu04_ a (LSCO) on the assumption that the 
electron-phonon coupling constant A ep has the same OIE 
on the in-plane effective supercarrier mass. Actually, 
Eq. (2) used in the Comment is incorrect. The correct 
T c formula in the (bi)polaron theory has the polaronic 
half bandwidth in front of the exponent rather than the 
phonon energy [3|. When the correct expression is ap- 
plied, the theory describes well the doping dependence 
of the isotope exponents in many cuprate superconduc- 
tors Q. Also, the tunneling experiments [3, @| have 
consistently shown that X ep is larger than 2.5 for op- 
timally doped Bi 2 Sr 2 CaCu 2 8+a (BSCCO). Therefore, 
the weak-coupling BCS-like formula does not hold in 
cuprates. A strong-coupling formula has been used to 
consistently explain the negligible OIE on T c and sub- 
stantial OIE on the in-plane effective supercarrier mass 
in nearly optimally doped BSCCO [fj|. In the under- 
doped regime, superconductivity should be better de- 



scribed by the Bose-Einstein condensation of local pairs 
in the strong coupling limit Q- In this case, T c is essen- 
tially proportional to n s /m*£ (where n s is the supercar- 
rier density and m*£ is the in-plane effective supercarrier 
mass), in agreement with the well-known Uemura plot 
Q. This implies that the OIE on T c is essentially pro- 
portional to the OIE on mil . This scenario can naturally 
explain why the magnitudes of the exponents for the OIE 
on both T c and m*j* increase with the decrease of dop- 
ing and are even larger than 0.5 in deeply underdoped 
samples (e.g., Lai.94Sr .06CuO4_ y ) [HQ. Therefore, the 
claim that the polaronic effects cannot explain the doping 
dependence of the OIE on T c is lack of scientific ground. 

Then the author of the Comment 0] attempts to ex- 
plain the doping dependence of the OIE on T c in terms 
of his own theory based on the t — J model. According to 
his model Q, AT C /T C = (1/A)AJ/J, where A = 0.2-0.3. 
Since the OIE on J was found to be —0.9% for undoped 
YBa 2 Cu 3 6 (Ref. 0), the predicted AT C /T C for the op- 
timally doped YBa 2 Cu 3 6 .94 should be -(2.7 - 4.5)%, 
in disagreement with the measured value of —0.27% 
(Ref. [ll|). In the deeply underdoped regime, T c oc I /mil 
and I /mil ^ J within the t — J model, so AT C /T C = 
AJ/J — —(0.6 — 0.9)%, in disagreement with the mea- 
sured values of -(5 - 12)% (Ref. [§]). Therefore, the 
t — J model cannot explain the doping dependence of 
the OIE on T c . Further, the author of the Comment Q 
argues that since the site-selective OIE shows that the 
OIE on T c mainly contributes from the planar oxygen, 
the apical oxygen is not important for the pairing. This 
argument is in parallel with the statement that since the 
optimally doped YBa 2 Cu30g.94 has a negligible OIE on 
T c , phonons are not important to the pairing mechanism. 
In fact, the site-selective OIE experiment [l2j shows that 
the apical oxygen contributes about 40% of the OIE on 
mH m YBa 2 Cu3 6.94- This suggests that the apical oxy- 
gen is important to the pairing, in agreement with our 
argument based on the bulk-sensitive x-ray-absorption 
experiment [l3| on (Yi_ 2: Ca a; )Ba 2 Cu307_ 2/ . 

The author of the Comment also criticizes the conclu- 
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sion about strong couplingto multiple phonon modes re- 
vealed by both tunneling [J, H, [I34l7l and angle- resolved 
photoemission spectra (ARPES) [17J, [l8[. There are sev- 



eral important facts about the strong coupling features 
revealed by both tunneling and ARPES data. First, the 
energies of strong coupling features match very well with 
those of the phonon modes revealed by the neutron 14 1 
and Raman [15J data. Second, the energies of strong cou- 
pling features revealed by tunneling spectra match very 
well with those of strong coupling features revealed by 
ARPES 17J. Third, the energies of strong coupling fea- 



tures in different cuprate systems such as LSCO, YBCO, 
and BSCCO are very similar IJ, llTJ, [l8| and nearly inde- 
pendent of doping [l9(. Such excellent consistencies un- 
ambiguously demonstrate that these strong-coupling fea- 
tures are intrinsic and arise from strong electron-phonon 
interactions. Furthermore, a detailed review of tunnel- 
ing, ARPES, and optical experiments has recently been 
given by Maksimov et al. These authors pro- 

vide consistent evidence for strong coupling to multiple 
phonon modes from the tunneling, ARPES, and optical 
results. 

Concerning the spin-fluctuation pairing, we have 
shown that the magnetic resonance mode revealed 
by neutron experiments plays a minor role in high- 
temperature superconductivity [H HH . One might argues 
that since the magnetic resonance mode is only a small 
fraction of the spin-excitation spectrum [20j j , its contribu- 
tion to the electron pairing should be insignificant while 
the coupling to the whole spin-excitation spectrum may 
still play an important role in the d-wave pairing. A re- 
cent theoretical calculation [2l| has shown that strong 
coupling to the whole spin-excitation spectrum of under- 
doped YBa2Cu306.6 (measured by neutron scattering) 
can lead to d-wave high-temperature superconductivity 
with T c = 174 K and the coupling constant = 1.39. 
It is important to note that the authors of Ref. [5l| have 
used a large renormalized coupling strength U (1.59 eV), 
which is too large compared with that estimated from 
several other independent experiments (see a recent re- 
view article These experiments consistently show 
that U < 0.17 eV. This implies that \ d < 0.0115 
(since Xd oc U 2 ). Such a small coupling constant cannot 
lead to high-temperature superconductivity. Further ex- 
perimental evidence for no magnetic pairing mechanism 
is that T c is very insensitive to the magnetic spectral 
weight, as clearly demonstrated from the neutron data 
[22I of slightly underdoped YBa2Cu30g.92 and slightly 
overdoped YBa2Cu306.97- The two compounds have al- 
most the same T c (91-93 K), but the magnetic spectral 
weight for YBa2Cu306.97 is at least three times smaller 
than that for YBa2Cu3 06.92- As pointed out by Maksi- 
mov et al. pj|, the little dependence of T c on the mag- 
netic spectral weight is incompatible with the magnetic 
pairing mechanism. On the theoretical ground, recent 
variational Monte Carlo simulations 23J, which are based 
on an advanced sign-problem-free Gaussian-Basis Monte 



Carlo algorithm, have shown that the simplest Hubbard 
model, advocated by Plakida and some other authors, 
does not account for high-temperature superconductiv- 
ity. 

Another important issue is the intrinsic pairing sym- 
metry in the bulk of superconducting cuprates. Because 
nearly all the surface and phase-sensitive experiments for 
both electron- and hole-doped cuprates provide clear ev- 
idence for d-wave order-parameter (OP) symmetry 24j . 
the d-wave pairing symmetry has become an indisputable 
fact to most researchers in the field. However, it is im- 
portant to note that these surface and phase-sensitive 
experiments based on planar Josephson tunneling are 
probing the OP symmetry at surfaces and interfaces, 
which were found to be underdoped 2^, 26 1. Since the 
majority of charge carriers are oxygen-hole bipolarons 
in the underdoped regime [l| and the OP symmetry of 
the Bose-Einstein_condensate of the oxygen-hole bipo- 
larons is d-wave [271 ] . the phase-sensitive experiments just 
probe the <i-wave OP symmetry of the dominant com- 
ponent. Since the OP symmetry of the Bose-Einstein 
condensate has nothing to do with the pairing symme- 
try, the phase-sensitive experiments do not probe the 
pairing symmetry associated with the pairing interac- 
tion. In order to probe the intrinsic pairing symmetry, 
bulk-sensitive data should be obtained from significantly 
overdoped samples where the dominant charge carriers 
are Fermi-liquid-like and the superconducting transition 
is mean- field- like [l[. Based on the quantitative analyses 
of many bulk-sensitive experiments (in addition to some 
bulk- and phase-sensitive experiments such as nonmag- 
netic pair-breaking effects), we have concluded that the 
intrinsic pairing symmetry in the bulk of cuprates is not 
d-wave, but extended s-wave (having eight line nodes) in 
hole-doped cuprates 1281 and nodeless s-wave in electron- 
doped cuprates [5^, |30| . 

The author of the Comment does not believe the intrin- 
sic pairing symmetry inferred from the bulk and phase- 
sensitive nonmagnetic pair-breaking effects. He argues 
against the extended s-wave gap symmetry using the 
surface-sensitive ARPES and Fourier transform scanning 
tunneling spectroscopy (FT-STM) experiments. Even 
surface-sensitive ARPES data of nearly optimally doped 
BSCCO can be better explained in terms of an extended 
s-wave gap 17,[2^,|3l|. The gap along the diagonal direc- 
tion Ad is small (< 7 meV) for nearly optimally doped 
BSCCO 1 171. 1 3 ill but becomes larger in heavily overdoped 
BSCCO [28|, |32j. Furthermore, both surface-sensitive 
ARPES and FT-STM data of a nearly optimally over- 
doped BSCCO can also be well explained in terms of an 
extended s-wave gap with ~ 4 meV (Ref. 31]). The 
significant uncerntainty in extracting the gap size from 
the ARPES data in a slightly overdoped and two under- 
doped BSCCO crystals [33j does not allow one to make 
distinction between a d-wave gap and an extended s-wave 
gap with Ad < 4 meV. In fact, no ARPES data along 
the diagonal direction were given in this study [33j j , which 
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FIG. 1: Photoemission spectrum along the diagonal direction 
for a nearly optimally doped Bi2Sr2Ca2Cu30io+i/ (Bi-2223). 
The spectrum was taken with an energy resolution of 10 meV 
(Ref. [I(J). The solid line is the curve calculated with the 
fitting method of Ref. [34| and the following parameters: Ad 
= 7.0 meV and V = 9.0 meV. 



makes it harder to draw a definitive conclusion about the 
gap symmetry from the ARPES data. 

In order to further prove our extended s-wave pair- 
ing symmetry in hole-doped cuprates, we determine the 
Ad value from some high-resolution ARPES data. The 
fitting method used in Ref. [HI for extracting the gap 
size from ARPES data should be the most reliable since 
the gap sizes extracted from this method match pre- 
cisely with those independently determined from the FT- 
STM data " 
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Fig. 1 shows photoemission spectrum 
along the diagonal direction for a nearly optimally doped 
Bi 2 Sr2Ca2Cu30io+ a (Bi-2223). The spectrum was taken 
with an energy resolution of 10 meV (Ref. [36]). The solid 
line is the curve calculated with the method of Ref. [33 ] 
and the following parameters: = 7.0 meV and V — 
9.0 meV (where T is the electron life-time broadening 
parameter). It is apparent that the gap size along the 
diagonal direction is not zero. 

In Fig. 2, we show photoemission spectra along the 
diagonal direction for a nearly optimally doped BSCCO 
(Fig. la) and a heavily overdoped BSCCO (Fig. lb). The 
spectra for the nearly optimally doped and heavily over- 
doped BSCCO crystals were taken with energy resolu- 
tions of about 6 meV and 10 meV, respectively (32l. 
The solid lines are the curves calculated with the fol- 
lowing parameters: Ad = 6.0 meV and T = 6.0 meV 
for the nearly optimally doped BSCCO; A D = 14 meV 
and r = 18 meV for the heavily overdoped BSCCO. It 
is striking that Ad increases from the increase of dop- 
ing, in agreement with the earlier ARPES result [35| and 
break-junction tunneling experiments [28| . 

In summary, our conclusions of the phonon-mediated 
pairing mechanism and s-wave pairing symmetry in 
cuprates are well supported by experiments. The crit- 
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FIG. 2: a) Photoemission spectrum along the diagonal direc- 
tion for a nearly optimally doped BSCCO. The spectrum was 
taken with an energy resolution of about 6 meV (Ref. [37]]). 
b) Photoemission spectrum along the diagonal direction for a 
heavily overdoped doped BSCCO with T c — 60 K. The spec- 
trum was taken with an energy resolution of about 10 meV 
(Ref. [1^). The solid lines are the curves calculated with the 
following parameters: Ad — 6.0 meV and T — 6.0 meV for 
the nearly optimally doped BSCCO; Ad = 14 meV and F = 
18 meV for the heavily overdoped BSCCO. 



icisms raised in the Comment Q are lack of scientific 
ground although the author of the Comment agrees that 
there exist polaronic charge carriers and electron-phonon 
coupling is strong in cuprates. 

* gzhao2@calstatela.edu 



[1] 



[2] 



[3] 
[4] 



Guo-meng Zhao, Pairing mechanism of high- 
temperature superconductivity: Experimental constraints, 
arXiv:1012.2368 , to be published in Physica Scripta. 
N. M. Plakida, Comment on the paper "Pairing mecha- 
nism of high-temperature superconductivity: Experimen- 
tal constraints", to be published in Physica Scripta. 
A. S. Alexandrov, Phys. Rev. B 46, 14932 (1992). 
R. S. Gonnelli, G. A. Ummarino, and V. A. Stepanov, 



4 



Physica C 275, 162 (1997). 
[5] D. Shimada, Y. Shiina, A. Mottate, Y. Ohyagi, and N. 

Tsuda, Phys. Rev. B 51, 16495 (1995). 
[6] G. M. Zhao, V. Kirtikar, and D. E. Morris, Phys. Rev. B 

63, 220506R (2001). 
[7] A. S. Alexandrov and N. F. Mott, Polarons and Bipo- 

larons, 67-106 (World Scientific, Singapore, 1995). 
[8] Y. J. Uemura et. al, Phys. Rev. Lett. 62, 2317 (1989). 
[9] G. M. Zhao, K. Conder, H. Keller, and K. A. Muller, J. 

Phys.: Condens. Matter, 10, 9055 (1998). 
[10] G. M. Zhao, Phys. Rev. B 75, 104511 (2007). 
[11] G. M. Zhao and D. E. Morris, Phys. Rev. B 51, 16487(R) 

(1995). 

[12] G. M. Zhao, J. W. Ager III, and D. E. Morris, Phys. Rev. 

B 54, 14982 (1996). 
[13] M. Merz, N. Niicker, P. Schweiss, S. Schuppler, C. T. 

Chen, V. Chakarian, J. Freeland, Y. U. Idzerda, M. 

Klaser, G. Muller- Vogt, and Th. Wolf, Phys. Rev. Lett. 

80, 5192 (1998). 
[14] G. M. Zhao, Phys. Rev. B 75, 214507 (2007). 
[15] H. Shim, P. Chaudhari, G. Logvenov and I. Bozovic, 

Phys. Rev. Lett. 101, 247004 (2008). 
[16] G. M. Zhao, Phys. Rev. Lett. 103, 236403 (2009). 
[17] G. M. Zhao, Phys. Rev. B 75, 140510(R) (2007). 
[18] X. J. Zhou et at, Phys. Rev. Lett. 95, 117001 (2005). 
[19] E. G. Maksimov, M. L. Kulic, and O. V. Dolgov, Adv. 

Condens. Matter Phys. 2010, 423725 (2010). 
[20] Hae- Young Kee, Steven A. Kivelson, and G. Aeppli, 

Phys. Rev. Lett. 88, 257002 (2002). 
[21] T. Dahm, V. Hinkov, S. V. Borisenko, A. A. Kordyuk, 

V. B. Zabolotnyy, J. Fink, B. Buchner, D. J. Scalapino, 

W. Hanke, and B. Keimer, Nature Physics 5, 780 (2009). 
[22] P. Bourges, 1998, The Gap Symmetry and Fluctuations 

in High Temperature Superconductors, NATO ASI Series, 

Physics, vol. 371, edited by J. Bok, G. Deutscher, D. 

Pavuna and S. A. Wolf (Plenum Press, New York, 1998), 

p. 349. 



[23] T. Aimi and M. Imada, J. Phys. Soc. Jpn. 76, 113708 
(2007). 

[24] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 
(2000). 

[25] J. Betouras and R. Joynt, Physica C 250, 256 (1995). 
[26] J. Mannhart and H. Hilgenkamp, Physica C 317-318, 
383 (1999). 

[27] A. S. Alexandrov, Physica C 305, 46 (1998). 

[28] G. M. Zhao, Phys. Rev. B 64, 024503 (2001). 

[29] G. M. Zhao, Phys. Rev. B 82, 012506 (2010). 

[30] G. M. Zhao and J. Wang, J. Phys.: Condens. Matter 22, 

352202 (2010). 
[31] G. M. Zhao, Phil. Mag. B 84, 3861 (2004). 
[32] I. Vobornik, R. Gatt, T. Schmauder, B. Frazer, R. J. 

Kelley, C. Kendziora, M. Grioni, M. Onellion, and G. 

Margaritondo, Physica C 317-318, 589 (1999). 
[33] W. S. Lee, I.M. Vishik, K. Tanaka, D.H. Lu, T. Sasagawa, 

N. Nagaosa, T.P. Devereaux, Z. Hussain, and Z.-X. Shen, 

Nature 450, 81 (2007). 
[34] H. Ding, J. C. Campuzano, A. F. Bellman, T. Yokoya, M. 

R. Norman, M. Randeria, T. Takahashi, H. Katayama- 

Yoshida, T. Mochiku, K. Kadowaki, and G. Jennings, 

Phys. Rev. Lett. 74, 2784 (1995). 
[35] R. J. Kelley, C. Quitmann, M. Onellion, H. Berger, 

P. Almeras, and G. Margaritondo, Science 271, 1255 

(1996). 

[36] D. L. Feng, A. Damascelli, K. M. Shen, N. Motoyama, 

D. H. Lu, H. Eisaki, K. Shimizu,3 , J.-I. Shimoyama, K. 

Kishio, N. Kaneko, M. Greven, G. D. Gu, X. J. Zhou, C. 

Kim, F. Ronning, N. P. Armitage, and Z.-X Shen, Phys. 

Rev. Lett. 88, 107001 (2002). 
[37] J. D. Koralek, J. F. Douglas, N. C. Plumb, Z. Sun, A. V. 

Fedorov, M. M. Murnane, H. C. Kapteyn, S. T. Cundiff, 

Y. Aiura, K. Oka, H. Eisaki, and D. S. Dessau, Phys. 

Rev. Lett. 96, 017005 (2006). 



